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The present invention relates to a method of determining filter coefficients 



from Line Spectral Frequencies comprising recomputing P(z) and Q(z) polynomials and 
comprising calculating the a) x coefficients. 



5 communications since the coded speech signal can be transmitted in a manner in which the 
redundancy commonly experienced in human speech is reduced. Linear Predictive Coding 
(LPC) is a known technique normally used in speech coding and in which the correlation of 
the speech signal is removed by means of a filter. The filter is best described by way of one 
of a different set of parameters, and one important set of which comprises LSFs. 
10 An accurate representation of the filter is an important requirement since such 

information is transmitted with the speech signal for subsequent reconstruction of the speech 
signal at a signal-receiving unit. 



have been well-documented since the inception of this concept in 1975. As is well known, 
15 the representation of an inverse LPC filter A(z) in the form of LSFs is derived from the 
representation of A(z) by its set of zeros in the z-plane. Insofar as the function A{z) 
represents an all-zero filter, it can be fully and accurately described by way of reference to its 
corresponding set of zeros. 



20 polynomial A m (z) of order m into two inverse polynomial functions P(z) and Q(z) . For 
confirmation, the polynomial A m (z) and the two inverse polynomials appear as follows: 



The coding of speech signals is used particularly in the field of mobile 



The advantages of representing LPC filter coefficients in the form of LSFs 



Computation of the LSFs commences with the decomposition of the 



A m (z) = l + a x z- x +a 2 z- 2 +„. + a m z 



-m 



25 



and 



Q(z) = A m {z)-z-^A m {z-') 
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The polynomials P(z) and Q{z) each have {m + 1) zeros and exhibit various 
important characteristics. In particular: 

5 

all zeros of P{z) and Q(z) are found on the unit circle in the z-plane; 

the zeros of P(z) and Q(z) are interlaced on the unit circle and the zeros do not overlap; and 

10 the minimum phase property of A m (z) is easily preserved when the zeros of P{z) and Q(z) 
are quantised. 

Analysis of the above confirms that z = -1 and z = +1 is always zero with the 
functions P(z) and Q(z) and since these zeros do not contain any information relating to the 

LPC filter, they can simply be removed from P(z) and Q(z) by dividing by (1 + z" 1 ) and 

15 (1-z- 1 ). 

Such revised functions can be represented when m is even as follows: 



_*2_ ande , (z)= _ew r 



20 and when m is odd as: 



The advantageous properties of functions P(z) and Q(z) as noted above are 
25 also valid for F{z) and Q\z) . Since the coefficients of F(z) and Q(z) comprise real 

numbers, the zeros form complex conjugate pairs such that the search for zeros only has to be 
conducted on the upper half of the unit circle, i.e. where 0 < co < n . 

It generally proves inconvenient to compute complex zeros, particularly by 
way of computerised numerical analysis methods, and so the functions P\z) and Q(z) are 
30 transformed to functions P"(z) and Q\z) with real zeros. Also, the functions F(z) and 
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Q (z) always have an even order and, since they are symmetrical, the functions can be re- 
written with real zeros to the following manner: 



F'(a>) = 2Z Pl "cos((m p -i)a>) 

«=0 

5 

gV) = 2i>,"cos(K-*>) 

1=0 

where p 0 "=l, p lx . Mf - l n = P ux .j*,-\ > p m /=^p m /,q Q n = l, 9 u ..^-i"=^i,2. m 9 -i\ 

q m ^ "= ~9m q 1 » where m p is equal to the number of zeros of P f (z) on the upper half of 

10 the unit circle and where m q is equal to the number of zeros of Q\z) on the upper half of the 
unit circle. 

When seeking the zeros of these functions, advantage can be taken from the 
form of the representations for P"(z) and Q"(z) due to the fact that the number of zeros to 
be located is already known. 

15 Importantly, and of particular relevance to the present invention, once the 

LSFs have been identified and employed as required, the recomputation of the LPC filter 
coefficients from the LSFs is required. While this stage represents a much less 
computationally intensive calculation than the computation of the LSFs from the filter 
coefficients as discussed above, problems and disadvantageous limitations are experienced. 

20 In particular, the values of intermediate coefficients can be disadvantageously high and this 
can lead to numerical problems even when employing floating point representations. 

Recomputing LPC filter coefficients a t from LSFs is much less 
computationally intensive than computing the LSFs from the filter coefficients. Each LSF 
co t , i — 0,1,..., m - 1 contributes to a quadratic factor of the form, 

25 1-2 cos(tf> f )z~ l + z" 2 . The polynomials P\z) and Q{z) are formed by multiplying these 
factors using the LSFs that come from the corresponding polynomial: 

F (z) = Yl (1 - 2 cos(fl> 2 , )z~ x + z~ 2 ) 

(=0 
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e , (2)=fl( 1 - 2cos ( a, 2, + ,)^ 1 +* -2 ) 

1=1 

The polynomials P(z) and Q(z) are computed by multiplying P\z) and Q(z) with the 
5 extra zeros at z = -1 and z = +1 . Finally, the filter coefficients are computed by using the 
following equation: 

10 which defines the relationship between the polynomial A m (z) and the two inverse 

polynomials discussed earlier. 

Thus, when recomputing the P{z) and g(z) polynomials one can use the 

above equations for P x {z) and Q\z) and add the extra zeros. Thus, for m is even: 

m p -\ 

15 P(z) = (l + z- l )Y[(l-2cos(0) 2l )z- 1 + z~ 2 ) 
Q(z) = (1 - z- l )fl(\ - 2cos(a> 2l+1 K 1 + z" 2 ) 

1=0 

and for m is odd: 

20 

P(z) = Y[ (1 - 2 cos(a> 2 , )z"' + z" 2 ) 

1=0 

= (1 - z- 1 )(1 + z^na - 2cos(o 2 , +1 )z- 1 + z" 2 ) 

25 Since <2> ( coefficients are ordered in increasing frequency the contributions of the first 

cos(*y, ) coefficients are positive and the last coefficients are negative. This introduces an 
undesirable growth of the intermediate coefficient values while conducting the polynomial 



^ * * «^HNL000385EPP 



5 03.07.2000 
multiplications (1-2 cos(fi> 2i+I )z" 1 + z~ 2 ) . With an increasing order of aw, such problems are 

amplified. To illustrate this, an example polynomial Q{z) = I — z~ 2N has been taken; it being 
noted that that m = 2N . Such a polynomial provides for a system with 2iV equidistant zeros 
on the unit circle. Since this is only a very simple example, it should of course be appreciated 
5 that in reality the growth of the intermediate coefficients can be much larger. Therefore a 
different strategy must be used. It has been found that for high order polynomials m=60 (or 
N=3Q) a double precision floating point representation will not be sufficient. A different 
method might involve searching for the best possible combination of co t and which exhibits 

the least amount of intermediate coefficient growth. However, due to the large number of 
10 possible combinations, this is most unlikely to be feasible and it also means that the optimal 
combination will never be found. 

The invention seeks to provide for a method of determining filter coefficients 
having advantages over known such methods. 

In accordance with the present invention there is provided a method of 
15 determining filter coefficients from Line Spectral Frequencies as noted above and 

characterised by the steps of addressing the polynomials in a series and reducing the number 
of polynomials in <o t in the said series by combining the polynomials in a> t two by two in a 

manner so as to arrive at two polynomials in co l and determining the product of the said two 
polynomials. 

20 The invention serves to combine a> t in such a way that hardly any signal 

growth occurs and proves particularly advantageous since the use of an increasing index i 
would not seem to offer a good solution. 

Using the method of the present invention with the example polynomial noted 
above, i.e. Q(z) = l — z~ 2N , the intermediate coefficients are never larger than 2. In practice 

25 only a limited amount of intermediate coefficient growth occurs. Advantageously, the 
invention need not comprise a particularly complex method. In general it requires only a 
different indexing and can advantageously deliver almost optimal results. For P(z) the same 
procedure could be used only if m is even then P(z) has a root at z = -1 . If m is odd P(z) 
does not have any additional roots, so no additional roots need to be added. 

30 Finally the relationship noted above is employed to compute the coefficients 

of A[z) from P{z) and Q(z) . 
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The invention is described further hereinafter, by way of example only, with 
reference to the accompanying drawing which is a graphical representation of the 
intermediate coefficient growth experienced in the prior art and in an example polynomial 
Q(z) = l-z- 2N . 



It is assumed that the original polynomial is reconstructed by combining the 
zeros with increasing co t . The maximum value of the largest coefficient during the 
10 recomputation procedure is plotted in the accompanying drawing. Note that the Y axis is 
logarithmic. For large order N the intermediate values of some of the coefficients become 
very high. 

However, such problems are not encountered in a method according to the 
present invention. 

15 As an example, and for Q{z) with m is even, the following ordering of the 

polynomials is used: 

v 0 [0] = l-z' 1 
v 0 [1] = 1 — 2 cos co x z~ l + z~ 2 
20 v 0 [2] = 1 - 2 cos <y 3 z -1 + z~ 2 

V 0 K 3 = 1 ~ 2 COS tfVm,-! 2 "' + Z ~ 2 

If m is odd then the terms for Q(z) are: 

25 

v 0 [0] = l-z" 1 

v 0 [l] = l — 2cosa>,z _1 +z~ 2 



30 



v 0 [m q ] = 1 - 2 cos a> 2 . mq _ x z' x + z 
v 0 [«, +l] = l + z _1 
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The next step is to combine the polynomials v 0 [i] . The strategy will be demonstrated with an 
example with m=12 and m q =6. The original seven polynomials are: 



vofOJ, v 0 [lj, v 0 [2], v 0 [3], v 0 [4], v 0 [5], and v 0 [6]. 

5 

In the first step the polynomials are combined two by two. Polynomial i is combined with 
polynomial [m q -i], this gives four intermediate polynomials v t [i]: 

v,[0] = v 0 [0]v 0 [6] 
10 v,[l] = v 0 [l]-v 0 [5] 
v,[2] = v 0 [2]-v 0 [4] 
v,[3] = v 0 [3] 

These four polynomials are combined in the same way, leading to two new polynomials 
15 v 2 fij: 

v 2 [0] = v 1 [0]-v 1 [3] 
v 2 [l] = v 1 [l]-v 1 [2] 

20 The product V2[0].V2[1] gives the final result: 
v 3 [0] = v 2 [0].v 2 [l] 

The procedure can formally be described by the following pseudo program: 

25 

if m is even 
begin 

niq—m/2 
m c —m q +l 

30 end 
else 
begin 
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m q =(m-l)/2 
m c —iriq+2 



end 



n p —m c 



5 i =mc> > 1 /*arithmetic shift right */ 
k=0 

while (i>0) 
begin 

n=0 

10 while (n<i) 

begin 

Vk+i[n]=Vk[nJ, v k [n p -n-l] 
n=n+l 

end 

15 if n p is odd then 

begin 

Vk+i[n]=v k [n] 
n=n+l 

end 

20 n p =n 
k=k+l 

i=n»L /*arithmetic shift right*/ 

end 



25 Using this method with the example polynomial Q(z) = 1 - z~ the 

intermediate coefficients are never larger than 2. In practice only a limited amount of 
intermediate coefficient growth occurs. This is not a very complex method (actually it only 
uses a different indexing) and delivers almost optimal results. For P(z) the same procedure 
will be used only if m is even then P(z) has a root at z = -1 . If m is odd P(z) does not have 

30 any additional roots, so no additional roots need to be added. The last step consists of using 
the equation 
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to compute the coefficients of A(z) from P(z) and Q(z) . 
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1. A method of determining filter coefficients from Line Spectral Frequencies 
comprising recomputing P(z) and Q(z) polynomials and comprising calculating the <d x 
coefficients, characterised by the steps of ordering the polynomials in a series and reducing 
the number of polynomials in w t in the said series by combining the polynomials in & t two 

5 by two in a manner so as to arrive at two polynomials in co l and determining the product of 
the said two polynomials. 

2. A method as defined in claim 1, wherein at least one series of intermediate 
polynomials is formed by combining the original polynomials two by two; the polynomials 

10 of the at least one intermediate series also being combined two by two so as to arrive at a yet 
further reduced number of polynomials. 

3. A method as defined in claim 1 or 2, wherein the following ordering of 
polynomials is used for m is even: 

15 

VotO^l-z" 1 

v 0 [l] = l — 2cos*»iZ' ! + z~ 2 

v 0 [2] = l-2cosfi> 3 z" 1 + z~ 2 

20 v 0 [m q ] = 1 - 2 cos co 2 *m q -i 2 1 + z ~ 2 

4. A method as defined in claim 1 or 2, wherein the following ordering of 
polynomials is used for m is odd: 

25 v 0 [0] = l-z- 1 
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v 0 [m q ] = 1 - 2 cos <y 2 . m<i _,z" 1 + z~ 2 
v 0 l>, + l] = l + z~' 



12 03.07.2000 

ABSTRACT: 



The present invention provides for a method of determining filter coefficients 
from Line Spectral Frequencies comprising recomputing P(z) and Q(z) polynomials and 
comprising calculating the co t coefficients, characterised by the steps of ordering the 
polynomials in a series and reducing the number of polynomials in to t in the said series by 
5 combining the polynomials in co t two by two in a manner so as to arrive at two polynomials 
in a> l and determining the product of the said two polynomials. 
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